A first estimation of chiral four-nucleon force effects in 4 He 
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Abstract 

We estimate four-nucleon force effects between different He wave functions by calculating the 
expectation values of four-nucleon potentials which were recently derived within the framework of 
chiral effective field theory. We find that the four-nucleon force is attractive for the wave functions 
with a totally symmetric momentum part. The additional binding energy provided by the long- 
ranged part of the four-nucleon force is of the order of a few hundred keV. 

PACS numbers: 21.45. +v,21.30.-x,25.10.+s 
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I. INTRODUCTION 

In a recent paper [j| the leading contribution to the four-nucleon force, V^n, has been 
derived within the framework of chiral effective field theory. It is governed by the exchange 
of pions and the lowest-order nucleon-nucleon contact interaction and includes effects due 
to the nonlinear pion-nucleon couplings and the pion self interactions constrained by the 
chiral symmetry of QCD. The individual pieces of corresponding to the diagrams in 
Fig. □ read Q 
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FIG. 1: The leading contributions to the four-nucleon force. Solid and dashed lines represent 
nucleons and pions, respectively. Graphs resulting by the interchange of the vertex ordering and/or 
nucleon lines are not shown. 
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Here, the subscripts refer to the nucleon labels and $ = Pi' — Pi with p^ and pi being the 
final and initial momenta of the nucleon i. Further, q 12 = <fi + qi = —q^ — q± = —^34 is the 
momentum transfer between the nucleon pairs 12 and 34. Diagrams (b), (d), (g), (h), (i), 
(j), (m), (o) and (p) lead to vanishing contributions to the four-nucleon (4N) force. The 
total short-range 4N force depends only on one low-energy constant Ct- 



II. CALCULATIONS 



We would like to estimate the magnitude of that 4N force in the 4N bound state. In order 
to simplify the calculations in a first attempt we assume that the momentum part of the 
4 He wave function is totally symmetric with respect to any permutations of the nucleons. 
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Thus we deal with the totally antisymmetric spin-isospin part | ^) of the total wave function 

I = ^= ({| S12 = Ml2 = >| S 3 4 = l,t 3 4 = >} S=0T=0 



{| S12 = 0,ti2 = 1 >| s 34 = 0,t 34 = 1 >} 5=0jT=0 ) , 



(2) 



where s^- and are the total two-nucleon subsystem spins and isospins. The curly brackets 
denote the coupling of the subsystems spins and isospins to the total spin (S = 0) and 
isospin (T = 0) of the 4N bound state. The state | £) can be expanded into the sum of 
product states 
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= ^j E s(i) I Xi(*)X2(i)Xs(i)X4(i)) | 771(1)772(1)773(^4(1)), (3) 

where Xj(i) iVji^)) i s the spin (isospin) state of the j th nucleon in the i th term of the sum, 
and s(i) denotes the sign of the i th product state. The "+" and "-" signs inside the kets 
stand for the +| and — \ spin and isospin projections, respectively. All the single nucleon 
states are normalized to 1 



(Xj(i) I Xj(i)) = I 
and consequently also the state | £) has the same norm 

(e 1 = i- 



(4) 



(5) 



The momentum part of the total wave function in the 4N center of mass (cm.) system 
depends on three relative (Jacobi) momenta 

P1-P2 

P = o 



2p 3 - (p x + p 2 
Q = 
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3^4 - (Pi +P2+ Pz) 



(6) 
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where Pi are the individual nucleoli momenta. Equations (JOj) can be inverted in order to 
express the individual momenta in terms of the relative momenta p, q and t: 

6p-3q-2t 



Pi 

P2 
P3 
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-6p- 3q- It 
6 

3q-t 



3 

p 4 = t (7) 

The assumption that the momentum part of the He wave function is totally symmetric 
is still very general and we make further restrictions. We assume that the momentum part 
can be written as a function of one variable, x, where 



x 



i & + p " + fi2 + = k ( p 2 + i f 2 + 1'" 2 ) • (8) 



which is the cm. kinetic energy of the 4N system, (m is the nucleon mass.) This implicitly 
means that we set all angular momenta to zero. We will later show to what extent this 
choice is realistic. Consequently we can write the full wave function | \&) as 

(pqt\^) = f(x)\0- (9) 

In order to calculate the matrix elements \ V^n \ ^) we calculate first the matrix 
elements in the spin-isospin space for the pieces of the 4N force given in Eq. ([TJ. For V a we 
consider first the following expression 

(£ I v i I = (Qi ■ Qui) (94 • 912) (£ I (#i ■ qi) (^4 • 94) ri • f 4 f 2 • r 3 I 

24 24 3 

= — (qi ■ m) (g 4 ■ 912) s{i)s(i)q 1 (a)q i ([3){xi{.3) I <? a | Xi(i)) 

»=1 j=l a,/3,7,(5=l 

(X2O') I Xz{i)) (Xs(j) I X2(«)) (XaU) I o-/3 I X4(i)) 
Mi) I ^7 I Mj) I r 5 I (773(7) I r.5 | 7/3(1)) (774O') | r 7 | jj 4 (i)) 

= (Qi ■ Q12) (?4 • m) (qi ■ qd ■ (10) 
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The intermediate multiple sums in Eq. (|10|) were obtained by means of the Mathematica 
program. In the same way we obtain the other expressions 

(f I V 2 I = - (Qi ■ Q12) (?4 • Q12) (f I (ffi • 9l) (<?4 • 94) n • f 3 f 2 • f 4 I 

= 3 (qi ■ q\ 2 ) (94 • ?ia) (?i • , ( n ) 

(f I I = (?i • 912) (f I (?i • gi) (dvt • <f 4 ) (912 x q 4 ) ■ a 3 (fj x f 2 ) • f 4 | f ) 

= -2 (ft ■ <f 12 ) (<f 12 x <f 4 ) ■ (ft x g 4 ) , (12) 
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(f I ^3 C I = (t I (*i ■ 9i) (^4 • 94) (912 x ft) • a 3 (fi x f 2 ) • f 4 10 

= 2 (ft 2 x ft) ■ (ft x ft) , (17) 

(f I I = (CI (^1 • 934) (?2 • 92) (?3 • 93) (^4 • 94) f • fa f 3 • r 4 I 
= (9*3 x ft) • (g 34 x ft) +2 (g 34 x ft) • (ft x ft) +5 (ft • ft) (g 34 • ft) , (18) 

(£ I V\ I = (£ I (<?i • 9i) (^2 • 92) (£3 • 9s) (^4 • 94) fi • f a f 3 • f 4 I 

= (93 x ft) ■ (ft x ft 4 ) + 2 (ft x ft) ■ (ft x ft) + 5 (ft • ft) (ft ■ ft) , (19) 

(f I Vl* I = (f I (^1 ■ 9i) ^2 ■ (91 x 9*12) (£3 x <r 4 ) ■ ft 2 fi • f 3 I 

= -2 (ft 2 x gl) 2 , (20) 
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(f I v 2 I = - {Qi • Q12) (f I (<?i • gi) (£3 x a 4 ) • gl 2 (n x f 2 ) • f 3 | f ) 

= 4 (gl • gl 2 ) 2 , 



(21) 



(£ I Vi I = it I (^1 ■ Qi) (^3 x £4) ■ 912 (fi x f 2 ) - fa |0 

= -4 (gi ■ q 12 ) , 



(22) 



(£ I I = (CI (*1 X ^2) • gi2 (<?3 X <7 4 ) • <?12 T 2 ■ T 3 I 

= -4g 12 2 - 



(23) 



Once the matrix elements ([10 p -([23j ) have been calculated we are left with the eighteen 
fold momentum space integrals. They can be written as 



I Van I 



dp fdq fdt fdp' fdq' [dt'g A4 (x)f(x)V i (p,q,t,p / ,q',t') f(x')g Ai (x'), (24) 



24 
(2^ 

where the functions V % (p,q,t,p',q',t'^ arise from introducing (JJJ) into (fTUj) - P3^) and the 
remaining expressions in (pQ) . The additional factors -7—^ and 24 arise from the wave function 
normalization and due to the fact that all nucleons' permutations yield the same result in 
the case of the totally antisymmetric wave function. The functions <?a 4 (x) and g Ai (x') are 
introduced since the expressions in (0) need to be regularized. We choose a simple form 

31 
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(25) 



so all the results will depend on the parameter A 4 . 

We consider two types of the 4 He wave functions. First one is a pure model Gaussian 

Q 

function |2| 

2 3/2 



III 



^9/4^9/4 eX P I p X 



(26) 



where the value of the parameter (3 is chosen after ,2] as 0.514 fm 2 . Wave functions of the 
second type are obtained in a quite different manner. We consider the wave functions which 

nn 

are solutions of the Schrodinger equation with the NLO chiral potentials [3j,|4j labeled by the 
following sets of the parameters (A, A): (400 MeV/c, 500 MeV/c), (550 MeV/c, 500 MeV/c), 
(550MeV/c,600MeV/c), (400 MeV/c, 700 MeV/c) and (550 MeV/c, 700 MeV/c). We 



wave function clq a\ a 2 

(400MeV/c,500MeV/c) 1.53266 40.4324 2.36626 12.5715 0.927233 
(550MeV/c,500MeV/c) 2.12619 86.6989 2.41787 14.4705 0.921551 

TABLE I: Parameters of the one dimensional fits 1)2 7jl for the two chiral NLO wave functions 
considered in this paper. 

checked that the wave functions with the same parameter A have very similar properties 
so the dependence on A is very weak. Thus we restricted ourselves to two cases only: 
(A,A) = (400MeV/c,500MeV/c) and (550 MeV/c, 500 MeV/c). For these two wave func- 
tions gained by rigorous solutions of the 4N Faddeev-Yakubovsky equations we extracted 
the component with the totally antisymmetric spin-isospin part. In both cases this compo- 
nent is dominant. It constitutes 94.3 % (88.7 %) of the original (400 MeV/c, 500 MeV/c) 
((550 MeV/c, 500 MeV/c)) wave function. Further we removed all contributions from the 
states with non-zero angular momenta. These components are small and represent only 0.3 
% and 2.9 % of the corresponding full wave functions. In this way we end up with the wave 
function components depending only on magnitudes of the momenta p, q and t, ^o(p,q,t), 
given on a certain grid. In order to facilitate the calculations, we represented *&o(p, q, t) by 
a one variable formula analogous to (|26[): 

f 2 (x) = {a + a lX a2 ) exp (-a 3 x a4 ) , (27) 

with the parameters a , ai, a 2 , a 3 and given in Table I. For the reader's orientation 
we show in Fig. El the components ^o(p,q,t) of the two chiral wave functions plotted as 
a function of x — ^ (p 2 + |g 2 + |t 2 ) together with the lines fitted according to (|2"Tj) . It 
is clear that the fits can be considered to be reasonable approximations to the underlying 
^o(p, Q, t) components only for small values of x. At larger x the values of ^o{p, q, t) are 
clearly underestimated. However, we assume in this first attempt that the main contributions 
to the expectation values come from the x region, for which the fits still reflect the bulk 
properties of the original ^o{p, Q, t)- That is why in the actual calculations we could use the 
simple analytical forms of (|27j) . Note that the very simple Gaussian wave function is close 
to the NLO fit with A= 400 MeV/c. 

In the practical calculations we used the basic Monte Carlo method and generated uniform 
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12 3 4 

x [fm" 1 ] 



FIG. 2: The values of *&o(p, q, t) for all possible combinations of p, q and t are plotted as a function 
of x = i (p 2 + |g 2 + |t 2 ^) for the (A,A)= (400 MeV/c, 500 MeV/c) case with black symbols and 
for the (A, A)= (550MeV/c, 500MeV/c) case with grey (cyan in color) points. The corresponding 
fits are represented by lines of the same color. The dashed line shows the Gaussian wave function 
(EH. 
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5(1) 
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MeV 
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ya 


-0.002906 


22 x 10" 


-6 


v c 


-0.005557 


25 x 10" 


-6 


v e 


-0.008462 


20 x 10" 


-6 


vf 


0.005692 


12 x 10" 


-6 


yk 


0.0005925 C T 


19 x 10~ 7 




V 1 


0.000622657 C T 


10 x 10~ 7 


Ct 


yn 


-0.000046044 C| 


55 x 10~ 9 


Crp 



TABLE II: Expectation values of the individual parts of the 4N force for the Gaussian wave 
function fi(x). The regulator function g\ 4 (x) with A4= 500 MeV/c is used. For the three last 
terms the value of the low energy constant (LEC) Ct in GeV -2 should be inserted. All the numbers 
should be additionally multiplied by the factor 24. 



distributions in each of the eighteen dimensions by means of the portable random number 

n 

generator ran2 from Ref. p. We found it sufficient to restrict the magnitudes of the relative 
momenta p, q and t to the following values p m ax= <?max= £max= 6 fm^ 1 . As primary tests of 
our Monte Carlo calculations we checked the norm and the internal kinetic energy of 4 He. 
These quantities can be calculated very precisely as three fold integrals but for tests were 
written as nine fold and (squared) even as eighteen fold integrals. The 4N force expectation 
values are approximated by 

r v N 

1 = jfdv^-Efi^ (28) 
i=i 

for which the one standard deviation error estimate reads 



6(1) = r 



N\ 



N / N 



E/ 2 W - E/W • ( 2 9) 

i=l \i=l / 

Here the points xi, X2, ■ ■ ■ , are uniformly distributed in the eighteen dimensional volume 
v. Tables II and III show our results for the Gaussian function fi(x) and the two first chiral 
NLO wave functions from Tab. I. We used N— 10 9 integral points. 

We show also in Fig. El the expectation values of the 4N force as a function of the Ct 
LEC. This is our final prediction, which includes all the required corrections. For Ct ~ 13 
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parts 


7(400, 500) 


5(1) (400, 500) 


7(550,500) 


5(1) (550, 500) 


of the 4N force 


MeV 


MeV 


MeV 


MeV 


ya 


-0.00434503 


17 x 10~ 5 


-0.00222788 
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v c 


-0.0084033 


19 x 10~ 5 


-0.00445691 


12 x 10~ 5 


v e 


-0.0133568 
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-0.00683624 
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vf 


0.00914028 
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0.00460722 


56 x 10~ 6 
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0.000926931 C T 


14 x 10~ 6 C T 


0.000489232 C T 


86 x 10~ 7 C T 


V 1 


0.000964454 C T 


77 x 10- 7 C T 


0.000512526 C T 


49 x 10~ 7 C T 


yn 


-0.00007243 Cf 


41 x 10~ 8 C| 


-0.00003812 C| 


26 x 10~ 8 Cf 



TABLE III: The same as in Tab. I for the two chiral NLO wave functions. Note that all the values 
should be additionally corrected for the norms of the wave functions: (\E' | ^>)= 1.093 ((A, A)= 
(400MeV/c,500MeV/c)) and | *)=1.011 ((A,A)= (550 MeV/c, 500 MeV/c)). As in Tab. I, the 
factor of 24 is not included. 



GeV 2 the magnitudes of the sum of the expectation values reach their minimum and we 
obtain approximately -0.077, -0.107 and -0.061 MeV for the three wave functions (Gaus- 
sian, (A,A)= (400MeV/c,500MeV/c), (A,A) = (550 MeV/c, 500 MeV/c)) considered, re- 
spectively. For Ct=0 the corresponding numbers are -0.270, -0.386 and -0.219 MeV. Only 
the two parametrizations of the chiral wave functions are consistent, at least to some extent, 
with the 4N potential. Thus we can state that the 4N force effects might vary from a few 
tens of keV to 1-2 MeV. Note that not the whole range of the Ct values shown in the figures 
actually appears for different orders of the chiral expansion "\. 

It remains to check the influence of different regulator functions on our predictions. To 
this aim we took the first chiral wave function and calculated the expectation values addi- 
tionally with A 4 = 200, 300, 400 and 600 MeV/c. As can be seen in Fig-IH the results do not 
differ much from each other for A4 > 300 MeV/c. Note that our definition of the regulator 
function^ (x) given in (J23j) introduces an additional factor of 2, as compared for example 
with Q,9]. Thus the values 300 and 400 MeV/c for A 4 multiplied by roughly correspond 
to the parameters A (400 - 550 MeV/c) of the wave functions. 
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-20 -10 10 20 30 40 
C T [1/(GeV) 2 ] 

FIG. 3: The expectation values of the 4N force for 500 MeV/c as a function of the Ct LEC 
for different parametrizations of the 4 He wave function. The dashed line represents the Gaussian 
wave function (|26|). the dotted line corresponds to the case of (A, A)= (550 MeV/c, 500 MeV/c) 
and the solid line is for (A,A)= (400MeV/c, 500MeV/c). 

III. SUMMARY 

We estimated for the first time 4N force effects in 4 He by calculating explicitly the expec- 
tations values of different 4N force parts between several 4 He wave functions. Our estimates 
agree qualitatively with modern nuclear force predictions for the a particle [t| , which do not 
leave much room for the action of 4N forces. Our predictions lack full consistency between 
the wave functions and the 4N potential and also neglect smaller components of 4 He. The 
strong dependence of the expectation value on Ct in the considered interval will probably 
be reduced using a fully consistent 4 He wave function at order NNNLO. Nevertheless, our 
results give some hint how important 4N force effects might be. 
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